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Abstract 

Deep space laser ranging missions like ASTROD I (Single-Spacecraft Astrodynam- 
ical Space Test of Relativity using Optical Devices) and ASTROD, together with 
astrometry missions like GAIA and LATOR will be able to test relativistic gravity 
to an unprecedented level of accuracy. More precisely, these missions will enable us 
to test relativistic gravity to 10" — 10~^ of the size of relativistic (post-Newtonian) 
effects, and will require 2nd post-Newtonian approximation of relevant theories of 
gravity. The first post-Newtonian approximation is valid to 10~^ and the second 
post-Newtonian approximation is valid to 10~^^ in terms of post-Newtonian effects 
in the solar system. The scalar-tensor theory is widely discussed and used in tests 
of relativistic gravity, especially after the interests in inflation models and in dark 
energy models. In the Lagrangian, intermediate-range gravity term has a similar 
form as cosmological term. Here we present the full second post-Newtonian approx- 
imation of the scalar-tensor theory including viable examples of intermediate-range 
gravity. We use Chandrasekhar's approach to derive the metric coefficients and the 
equation of the hydrodynamics governing a perfect fluid in the 2nd post-Newtonian 
approximation in scalar-tensor theory; all terms inclusive of 0(c~^) are retained 
consistently in the equations of motion. 
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1 Introduction 



Although Einstein's general relativity has achieved great success both in ex- 
perimental tests and in astrophysical applications during the last few decades, 
the desire to find a gravitation theory consistent with quantum theory together 
with the ever-increasing precision of experiments and astrophysical observa- 
tions has urged many "alternative theories" to be proposed. Among them, the 
scalar-tensor theory is the most eminent one, because it is the simplest and 
most natural way to modify general relativity. Many modern theories, such as 
extra-dimensional theory, string theory, brane world and noncommutative ge- 
ometry, which try to unify gravity and microscopic physics or explain the dark 
energy in cosmology, demand a scalar field in addition to the metric tensor. 
In this paper, we take a phenomenological point of view. In the low energy 
effective field limit, the action of most of these theories can be transformed to 
the following general form: 

5 = ^/(^0i?-^0'>,. + 20A(0)-^£,(^7/..,0,^))v^d^^, (1) 

where 6{(f)) and A(0) are two arbitrary functions of the scalar field (p. g = 
det (gf^u) < is the determinant of the metric tensor (yf^,^, R is the Ricci scalar, 
\1/ denotes all the matter fields. Greek indexes run from 0, to 3, and Latin 
indexes run from 1 to 3. The signature of g^i, is (— , +, +, +). 

In general, the matter fields \l/ interact with both the metric field and the scalar 
field. If we assume the Einstein equivalence principle to be correct, then the 
matter fields \1/ do not interact directly with the scalar field and the interac- 
tion Lagrangian Ci{g^u, (p, ^) becomes Ci[g^y, '^). Since Einstein equivalence 
principle is verified to a very high accuracy (Ni, 2005), we shall assume it is 
valid here. Violations of the Einstein equivalence principle have been consid- 
ered in Ni (2005). Here, for simplicity, we do not consider it. Therefore, in this 
paper, we adopt the following action 

s = ^j[^R- ^-J^"^'^"^'- + - ^^lia,., ^)) V^d'x. (2) 

Depending on the functional form, the A(0) term can include (i) cosmologi- 
cal constant and quitessence; (ii) the mass term of scalar field which induces 
intermediate-range gravity. Since this paper deals with post-Newtonian ap- 
proximation and assume an asymptotic fiat spacetime, we include the case 
with the mass term and intermediate-range gravity only; we assume that 

A(0) = A2(0-0o)', (3) 
where A2 is a constant and 0o is the constant asymptotic value of (p. In this in- 
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vestigation, we will obtain the second post-Newtonian approximation of scalar- 
tensor theory of gravity including intermediate-range gravity using Eq. (2) 
with A(0) given by Eq. (3). 

Fujii and Maeda (2005) on the scalar-tensor theory of gravitation gives a 
good account of the historical development of the scalar-tensor theories. Here 
we present a very brief history related to our choice of Lagrangian/action 
Eq. (2). Jordan (1955) first proposed scalar-tensor theory in connection with 
projective geometry and five-dimensional Kaluza-Klein theory (Kaluza, 1920; 
Klein, 1926). Assuming the validity of Einstein equivalence principle, Brans & 
Dicke (1961) reached a specification of Jordan's theory. This theory is termed 
Brans-Dicke- Jordan theory in the compendium of metric theories of gravity 
compiled in Ni (1972). This theory is a special case of (2) with ^(0) = uj = 
const, and A = 0. Bergmann (1968) generalized Jordan's theory regarding to 
0. In his paper, the interaction Lagrangian included only electromagnetic field 
F^y, not other matters, with 



C-i = /(0)M, (4) 

where M is the Maxwell scalar formed from F^^, and / is an arbitrary function. 
From experiments on the test of weak equivalence principle, /(</>) is constrained 
as follows (Ni, 2005): 



^ < 10-^", (5) 

where U is the Newtonian potential. With /(0) = 1, Bergmann's theory is of 
the form (2) with denotes electromagnetic field F^^^. Wagoner (1970) com- 
pleted Bergmann's theory with other matters and postulated his principle of 
mutual coupling (equivalent to Einstein equivalence principle). The resulting 
theory is called Bergmann- Wagoner theory in the compendium of metric theo- 
ries of gravity compiled in Ni (1972). Bergmann- Wagoner theory is equivalent 
to Eq. (2). 

Deep space laser ranging missions such as ASTROD I (Single-Spacecraft As- 
trodynamical Space Test of Relativity using Optical Devices) and ASTROD 
(Astro dynamics Space Test of Relativity) (Bec-Borsenberger et al., 2000; Ni, 
2002; Ni et al., 2004), together with astrometry missions such as Global As- 
trometric Interferometer for Astrophysics (http:/ /www.esa.int/esaSC/1 20377 
_index_0_m.htm) and Laser Astrometric Test of Relativity (LATOR) (Tury- 
shev et al., 2004) will be able to test relativistic gravity to an unprecedented 
level of accuracy in the solar system. More precisely, these missions will en- 
able us to test relativistic gravity to 10"^ — 10"^ of the size of relativistic 
(post-Newtonian) effects, and will require 2nd post-Newtonian approxima- 
tion of relevant theories of gravity. The first post-Newtonian approximation 
is valid to 10~^ and the second post-Newtonian is valid to 10~^^ in terms of 
post-Newtonian effects in the solar system. 
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Here we present the full second post-Newtonian approximation of the scalar- 
tensor theory (2) with A(0) given by Eq. (3) treating A2 as an independent 
parameter. We derive the metric coefficients (in Sec. 2); when A2 = 0, our 
result agrees with the result of Damour & Esposito-Farese (1992, 1996) with 
one scalar field where they are comparable. Damour and Esposito-Farese use a 
field-theoretical point-particle approach; we use Chandrasekhar's perfect fluid 
approach. In Sec. 3, we derive the equation of the hydrodynamics governing 
a perfect fluid in the 2nd post-Newtonian approximation; all terms inclusive 
of 0{c~^) are retained consistently in the equation of motion. In Sec. 4, the 
various conserved quantities to 0{c~^) are isolated with the aid of the energy- 
momentum complex. In Sec. 5, we discuss solar-system dynamics and test of 
scalar-tensor theory. In Sec. 6, we present an outlook for further works. 



2 Second post-Newtonian approximation 



2.1 Field equations 



Variation of the action Eq. (2) with respect to gap yields 



where □g(-) = {■)-af3g'^^- T^y is the stress-energy- momentum tensor of matter 
deflned as (Landau & Lifshitz, 1962) 

2^ ^ dg^^- dx'^ dg%' ' ^ ' 

and T is the trace of T'^'^. Here we write T^'^ in the following form 

c'T^"' = pic" + n)u''u'' + 7r'^^ (8) 

where p and H are the density and the internal energy in the co- moving frame, 
ti'^ is the dimensionless 4- velocity and vr^*^ is the anisotropic stress tensor. For 
perfect fluids, 

where p is an isotropic pressure. 

Variation of the action with respect to yields 



In the following, we assume Eq. (3), that is, A(0) = A2(0 — 0o) 
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2.2 Perturbation of the scalar field and the metric 



We assume that the scalar field can be expanded in power series around its 
background value 0o ^-s Kopeikin & Vlasov (2004) and define 

= 0O(1 + C), (11) 

where ( is dimensionless perturbation of the scalar field around <^o- 

In particular, decomposition of the coupling function 9{(f)) can be written as 



e{(j))=Uo + UJ,C+^U2C^ + ..., (12) 



where Uq = ^(0o) and cu„ = (d^^/dC^^) 



(j>=4>o- 



Following Chandrasekhar (1965, 1969b), we look for solutions of the field equa- 
tions in the form of a Taylor expansion of the metric tensor and the scalar 
field with respect to the parameter e = 1/c such that 



(2) (4) (6) 

goo = -l+e^hoo + e^hoo + €^hoo^ , (13) 

(3) (5) 

goi = e^hoi + e^hoi + -- - , (14) 

.(2) ,(4) 

9ij = Sij +e hij + e hij ^ , (15) 

(2) (4) 

( = e\ +... , (16) 

and 



(0) o (2) , (4) 

Too = Too + ^'Too + e'Too + ■ ■ • , (17) 
(1) 3(3) 



Tni — eToi + e^Toi + ■ ■ ■ , (1^ 



(2) (4) 

Tij=e% + e% + --- . (19) 



Furthermore, we simplify the notations with the definitions: 

_ (2) _ (4) _ (3) _ (2) _ (2) 

N = hoo, L = hoQ Li = hoi Hy = hij, H = hkk, (20) 

_ (6) _ (5) _ (4) 

Q = hoQ, Qi = hoi, Qij = hij. (21) 

According to Eqs. (11), (12) and (16), we have 
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(2) /l (2)^ (4)\ 



(22) 
(23) 



0(e") means of order e'^. From here on, we omit the writing of O (£:"■) where 
there is no ambiguity. We introduce three parameters 7, /5 and i as follows: 



7 



a;o + 1 



cuo + 2 



(2cJo + 3)(2cJo + 4)^ 



_i (7-ir 

^ = Uo. 

2 7 + 1 



(24) 
(25) 

(26) 



For scalar-tensor theory with A2 = 0, it turns out that 7 and f3 become the val- 
ues of corresponding standard PPN parameters 7 and (3 in Will & Nordtvedt 
(1972). 

Corresponding expansion of the energy-momentum tensor is 



00 = p +e p 



-v^ + U- N --H 
2 2 



I 4 * 
+ e p 



1 

—'1 
4 



-v^ + -v'^U - -v'^H + v'^— - nA^ - -UH - -NH 



p^ 



1. 
2' 



1 

2' 



+ ^HikHik + ]^Hikv\^^ -L- LkV^ - ]^Qkk 



Toi = -ep*v' - e^p* 



V2 2 p* 



, (27) 

(28) 
(29) 



Here we have used the invariant density p* = ^/^vP p (Fock, 1959) with its 
2PN expression 



P* = P + e^p\^v^ + \h 



+ e^p 



8 2 4 2 

— -;HikHik + -H + LkV + -Q^fc 
4 o z 



(30) 
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2.3 The gauge condition 



We use the gauge condition imposed on the component of the metric tensor 
proposed by Kopeikin & Vlasov (2004) as follows: 




Although it is called "Nutku gauge" by Kopeikin and Vlasov, the gauge con- 
dition is different from Nutku (1969a,b). To 2PN order, this gauge gives 



1 (2) 

-H^i — -N^i — Hik^k + C 



1 1 1 1(2) 

— -^HikHik^i — -L^i + Li^t + -^Qkk,i — Qik,k — 2 C 



1(2) 



(2) (2) 



(4) 



+ 2 C - Hik C,fe - C Hik,k + C , 



(32) 



and 



,/ 1 1 (2) \ 

^'[-^N,--H, + Lk,k- Ctj 

+e' - ^NH^t - NN^t + \HikHik,t - HikU,k + NLk,k + \LkN ^k 
1 1 1 

+ -LkHk — LiHik^k — -L^t + Qk,k — -Qkk,t 



(2) 



(2) (4) 



1 (2) (2) 1 (2) _ _ _ 

-CNt-NC^,--CHt+ CLk,k + LkC,k- C,J =0. (33) 



2.4 Metric coefficients 



With the metric, Eqs. (13)-(15), the energy-momentum tensor, Eqs. (27)-(29), 
the decomposition of scalar field, Eqs. (11), (16), (22) and (23), the parameters, 
Eqs. (24)-(26) and the field equations (6) and (10), the metric coefficients can 
be solved. 



7 



2.4-1 Newtonian approximation 
The equation for N is 

(2) 

V^A^ = -SttGp* - 6 C, (34) 

(2) 

where ( satisfies 

J2) (2) 

V^C +eiC =-4(l-7)7rGp*, (35) 



and 



(2) ^ 

The solution for C (r ) is 



ei = 4^A20^. (36) 
1 + 7 



C in = (1 - 7)G / -/I (37) 

Hence, the solution for N{f) is 

(2) 

iV(r ) = 2G f P^df' + ^ ( ^df'. (38) 

(2) 

When A = 0, C ^ind reduces to (1 — 'j)U and 2f/, where 

U{f) = f -C^df', (39) 
Jv r — r 



is the Newtonian potential. 



2.4-2 1st post-Newtonian approximation 

Following the method developed by Chandrasekhar (1965, 1969b) and Nutku 
(1969a), we work out the 1st post-Newtonian approximation: 

(2) 

V^i/., = (5,,(-877rGp* + 6C), (40) 

and 

V^L, = 8(7 + l)7iGp*v\ (41) 

where G = 2/[0o(l + 7)]- According to Eq. (40), we have Hij = V5ij- For L, 
we have 



V^L = -87rGp* 



/ / p* 7 — 1 
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where 



1 1 (2) Aid - 1) (2) (2) 

(2) (2)2 (4) 



^3 = 4 — ^ A20O- 



7^ 



(4) 

The Poisson equation for ( is 



(4) 



l(7-iy + ^(7-i)v^-(7-i)n 

+ 3(7-l)4 + ^f 
p* 7 — 1 

TV V ? 



(7 - 1)^ 



with 



(2) (2)^ 



6 



3^2 _ _ IQ^ + 19 
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2-1 



From these Poisson equations, the integrals can be readily written out. 
2.4-3 2nd post-Newtonian approximation 

As above, we solve for the 2nd post-Newtonian approximation: 



(42) 
(43) 



(44) 
(45) 



-8(l+7)7rG'p*t;V 

-^-N,N, + VN,,-VV,,-W,V, 

(2) 1(2) 1(2) 1(2) 1(2) 

- C N,, - - C - -QN, + - C + 2 C 

(2) (2) (2)(2) 2{2-i - 1)(2) (2) 



+5. 



7-1 

11 D 

+ -7i;2 ^ _ + (27 - 1) — 



+ 



2 ' 2 

4/3-4 + 72-7(2) 



7-1 



C 



1 1 (2) 4f/3-l)(2) (2) 
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(2) (2)2 (4)^ 



(46) 



(7 + ^y{l^'' +^+ly+^-^(]+L^ 



11 1 
— -VikLk + -VkLi^k + ViLk,k — VkLk,i + Li^tt 



(2) 



1(2) 



1(2) 



(2) 



(2) 



3(2) 



1(2) (2) (2) (2) (2) 

t; C ,t^,i - C ,iLk,k - C Lk,ki - Li,k C ,k- C ,ikLk 



2(2-f - 1)(2) (2) (2)(2) 



7-1 



(2) 



(47) 



and 



V2Q = -87rGp*\ i(47 + 3^ + ^(27 + lyV + ^(27 + l)v2n 
18 4 2 

+ (7 + l)v^— -UN- -UV - 3-fN— + 3-fV— 
p* 2 p* p* 

+ lNV + lv''-L-Lkv''-l-Qkk 



+ 



2(7-1) 



(7 - 1)^ 



7-1 



4(/3-l)(7^ + 8/3-27-7) 
(7-1)3 



(2r 
C 



_ 3(72 - 7 + 4/3 - 4) (2) p ^ 4/5 - 7 - 3^1 
7-1 p* 7-1 J 

+2iViV,„ - ^iViV,fciV,fc + 2NtN,t - ^VNkVk + VN^u 
+3NVu + ^NtVt + VNkNk - NkL,k - Lk,kNt - 2NLk,kt 
-N,kLk,t - LkN^kt + N^kQki,i + N^kiQki - ^NiQkk,i - ^VtVt 
+L^u — ■^V,kL^k + '^VtLk,k — Vkt^k ~ Lk,iLi^k + Li kLi^k 
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(2) (2) (2) (2) (2) (2) 

(2) (2) (2) (2) (2) 

+2 C ,kLk,t — 2 C Lk^kt — 2 C ,t-^fe,fc — 2 C ^kt^k — 2 C ,fc-Z^,fct 
{2)(2) ^ ^(2)(2) ^ 2(27^-37 + 2/3- 1) (2) (2) 



+ C C,fciV,fc + 2C Ctt + ^^^r^TT^^^^ -C,tC,t 



_ 1) (2) (2) 2l (2)(2) (2) (4) 



4(/3-l)(7^ + 8/3-27-7) (2)(2) (2) 8(/3-l) (2) W 

(^^ — 1)"^ (7 — 1)2 ^ 



(2) (2)(4) (2)2 



+UVN + L) c - 2^3 c c + e3(iv - y) c 

(2)3 (4) (6) 



(6) 

where ( satisfies 



(6) 



7-1 4 

1; 



2_J^2n + 2_i,,2^2iv^^) 



4(/3 - ^ 7. 



7 



7 - 1 



p* 7 — 1 p* 8 



4(/9 - 1) ^(J) 



7-1 
2 



2(32/3^ + 67 -647 -6 + 32) (2)" 8(/j- 1)W| 



(7 - 1)= 



C + 



7 



c 



J 



(2) 1(2) 3(2) 1 (2) 1 (2) 



(2) 1 (2) (2) 



2 

(2) (2) 

C ,kLk,t + C ,kQki,i 



(2) 1 (2) 



(2) 4(/3-l) (2) (2) 

^ ^ + _ 1)2 ^ ,t 



32(^-l)2 + 2^(7-l)(2)(2) (2) 8(/?-l) (2) (4) 



(7-1)^ 



(7 - 1)^ 



1(4) 1(4) (4) {2f {2f 

+^C,kNk--C,kV,k+ C,« + 6^C +e4C 

(2) (4) (4) 

+e5C c -ei^c, 



(49) 



with 
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Table 1 

Summary of the parameters. 



Parameter 


What it measures, relative to GR 


Value in GR 


Value in STT 


7 


How much space curvature {gij) 
is produced by unit rest mass? 
(Misner et al., 1973) 

How much 2nd order nonlinearity is there 
in the superposition law for gravity (ffoo)? 
(Misner et al., 1973) 


1 
1 


1^0+2 

1 1 1^1 

(2tJo+3)(2aJo+4)^ 


i 


How much 3rd order nonlinearity is there 
in the superposition law for gravity (s'oo)? 







2(3+2a;o)(c^O+2)3 



(7-1)3(7 + 1) 



^5 



^4 _ 4^3 _ 24^^2 ^ 30^2 ^ ^28/3^ + 48/37 
+ 467 - 527 - 4t - 280/5 + 153 j , 
4(37^ - 67 - 16/3 + 19) 



r 



-A2</>, 



0' 



and 



^6 



K2/372 - 272 - 32/32 - 4/37 - i7 + i + 47 + 66/3 - 34) 



(50) 
(51) 
(52) 



(7 - 1)3(7 + 1) 
2.4-4 Summary of the parameters 

The physical meaning of A2 is that it gives an inverse range ^1 (= 4(1 — 
7)A20o/(l + 7)) of intermediate gravity. With A2 7^ the theory violates 
the inverse square law for gravitation (Fischbach & Talmadge, 1998; Li & 
Zhao, 2005). When A2 = 0, 7 and /3 reduce to the standard PPN parame- 
ters in 1st post-Newtonian approximation. It is clear that, besides them, only 
one parameter t emerges in the 2nd post-Newtonian approximation, which 
represents the 3rd order nonlinearity in Qqq. Table 1 gives a summary of the 
parameters involved and their values in general relativity. 

In 1992, Damour and Esposito-Farese proposed a multiscalar-tensor theory 
(Damour & Esposito-Farese, 1992). When only one scalar field involved, the 
action reads, 

= ^g^^ J d'^Xy/g^[R^ - 2gi^''d^(pdu^p] + Sm['4'm,g^,u] (53) 
where and (7*^ are the physical metric and Einstein-frame metric respec- 
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Table 2 

A parameters' comparison between Damour and Esposito-Farese's (DE's) and ours 
when only one scalar field involved, where 7 = 7— 1, /? = /?— 1, ao = dlnA{(p())/d(po, 
Po = da{ifQ)/d(po, (3q = dl3{ipo)/difQ (see Damour & Esposito-Farese (1996) for 
details) . 



Order (PN) 


Paramters 


DE's 


Ours 


1 
1 

2 
2 


7 
Cde 




1 


1 l^ootl 


1^0+2 


2 (l+ag)2 


4(2LJo+3)(a;o+2)2 


(1+^2)3 


7 

_M! 




7 



tively, and (7^,^ = A^{Lp)g*^^. Hereafter, tilde will be dropped for clarity. After 
setting 

and 

3 1/dlnA' 

0. In a follow- 



(55) 



2\ d<p ) ' 

the action (53) could reduce to ours (2) in the case that A2 
ing paper (Damour & Esposito-Farese, 1996), Damour and Esposito-Farese 
derived 2PN approximation of their theory, and showed that it would intro- 
duce only two new 2PN parameters Sde and C,de- (Here we use the subscript 
or superscript "DE" to denote Damour and Esposito-Farese's results.) When 
there is only one scalar-field involved, Cde-, depending only on 7 = 7 — 1 and 
/3 = /9 — 1, is not a new parameter in 2PN; is the only new independent 
parameter in 2PN related to our parameter l by Ede = —l — 12/5^/7 (see Table 
2). Damour and Esposito-Farese derived the 2PN deviation from GR of the 
metric component Qqq, but not other metric components. Their 2PN deviation 



from GR SgBf is given by 



00 



^_DE 2 s 

3c6 



£de 



47iGaU^ + 



7 ^' 



01 



(56) 



Our corresponding equation is 



00 



With the help of 
and 



(57) 

(58) 
(59) 

Eqs. (56) and (57) can be transformed into each other and agree. This is a 
consistency check for own calculation. 



1. 



V2f/3 = UV^U'^ - U'^V^U. 



2rr2 tt2^2i 
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3 Equations of Motion 



The equations of motion for perfect fluid are derived in 2nd post-Newtonian 
approximation. 

T^y = yields the momentum equation, 



V ■ V 



2 p* 

,2 V 



+ -v^Ii + v^— + v^\N --V 



P \ 

N-^-v]+^N+'-^V 
2 J p* 8 



4 



_ ^NV + N^ + L + Lkv^ - ]^Quk 



+P,i - ^pN,i 



-nP-v^ + Iu+'-^ + In-'-v] 
' V4 2 2p* 2 4 ; 

P 2 



2 • 



+£^pl -v'—V v'—N-\ Li Li-Lkiv" 

^\ 2 dt 2 dt dt 2 ' 

-^N,{N + v') + ^V,{N-v')\ 



+e^p*< 



U + 2N - -V \ —N 
2 dt 



+ -\v' + U + N--V\—V + -v'{ —N - 5—V 



dt 



dt 



dt 



nJ ly^ + ^y^u + -v'N--v'V + -UN 
''V 16 4 2 8 2 

_ ^UV - -^V + -N' - -NV + -V' 
4 2p* 2 4 16 
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+ -L + -Lfct;^--gfcfc)] 

d 111 



;ives the continuity equation, 



V ■ V 



p* + e'^p* 



\y^ + li + N--V 
2 2 



1 4 * 

+ e p 



8 2 V 4 / p* 



V ■ [pv) + p* 



+ L + 2Lfet;'= - -Qfcfc 
2dt dt 2 ' y 



V ■ {pNv) + LkP,k 



-v^ + Yi- -V + 2N\—N 

2 2 Jdt 

3 2 3^ 3,, 15, \ d 

-2 -n + -iv — F ' 



4 



4 



dt 



2p* 4 , 



2 p* 2 



0. 
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4 Conservation Laws 



We adopt the energy-momentum complex to obtain tlie conservation laws in a 
simple way and we set A2 = for simplicity. The energy-momentum complex 
is defined as in Nutku (1969b) 



where 



in which 



fj,auf3 



2(-^7)0 



2^, a/3 



Xi^" = Ri^'' ---n^"'R- 



and 



(62) 
(63) 

■^^^ - g^''U,<py (64) 
(65) 



Kopeikin-Nutku gauge Eq. (31) gives 

f/,r/ = <PV^g''^<PV^ggn,<.p - {<PV^gg'^)A<PV^99n,P- (66) 

To obtain the order we need, we extend the metric expansion of g^u as 



goi = £^Li + e^Qi, 

9ij = Sij + 2e'^'j5ijU + e^Qij + e^Sij, 
with Sij defined by Eq. (69). 
Thus, we can obtain conserved quantities: 
(1) energy: 



(67) 
(68) 
(69) 



(2) linear momentum: 

(3) angular momentum: 



E = c2|(e°°-pn°v^)d3a;; 



(70) 
(71) 
(72) 



Modulo divergence (Chandrasekhar, 1969a), and can be expressed as 
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I 4 * 

+e p 



^i;^ + f 37 + ^ j v^U + ^t;^n + (47 - l)IiU 



^(2672 + II7 + 32/5-37 



32 

7 + 1 



2 p p* 



+e- 



4nG 



2 37^-107-3/5 + 7 



+ ^UtUt - ^lUUu 



16 



7 + 1 



\UU,U 



(73) 



— TTj + 



(74) 



p*v' + eV*^' f^^'' + n - f/ + 



A-kG 



UiUt 



+e'p*< 



^v'' + \-i + ^WU + ^v'li - ^(27 + ifU'^ 



-f/n + (37-l) — f/ + — i;' 

p* p* 

1 t 1 

- -L + Lfct;'^ + -Qkk + C 

I9 1 — 7^ 27 — In 

2 



2(1 + 7) 7 + 1 p* 



97 



4(1-/3) 8/3 ^ 
1-7 7+1 



If/ 



+£^ 



1 



SttG 



2 117^-67 + 4/3 + 11 



16 



7 + 1 



\UU,Ut + ^-UtL,i 
7 + 1 



2 4(7-1) 7-3 



7 + 1 



7 + 1 



7 + 1 



7 + 5 



-UL.^tt + (67 - 3)L,f/,fcf/, 



7 + 

4(27' + 7-2) 
7 + 1 



UUhLk, 



7 + 1 
2(/3-l)(7-3) 



2 1 

LihLk,t H —rLkL ki 



1-72 
1 



ULkUk^ 



7 + 1 



7 + 1 



7 + 1 
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7-1 2(7-1) 37-1 

-U^kQk,i — — U^iQk,k —-UQik,kt 



7+1 ' ' 7+1 

27 2 

— —7Qik,tUk H --rQikU kt 

7+1 7+1 

7-1 1 

— —QkkU^it H — 

7+1 7+1 

1 



7 + 1 
2(27-1 



7 + 
1 



ik,k 



7 + 1 



SkkM + 



7 + 1 



kM 



7 + 1 

2(7-1) (4) 



LiQik, 



ki 



7 + 1 



(75) 



and 



47rG 



2(7 + l)5ikUU^t + - UL 



k,i 



+8" 



SttG 



47(372-47 + 1) , 157^-147 + 3 



7 + 1 
(37-1)2 



Hk 



2(7 + 1) 



i,k 



(7 + I) 



7 + 1 



2(7 + 1 



H —{LiQiki — QiiLki + LiQki^i — LnQki + LiQuk) 

7 + 1 



37-1 
7 + 1 
1 



^ + ^^^-^5rk{u^tQii - UQu,t) 



7 + 1 



7 + 1 



2(7 + 1) 



2(7 - 1) 47 
-^-^'.fc H — ; — QiU^k H —rUQi^k 



7 + 1 



7 + 1 



(76) 



5 Solar System Dynamics and Test of Scalar-Tensor Theory 



To test the scalar-tensor theory including intermediate-range gravitational 
force in the solar system, we need equations of motion for solar-system dy- 
namics. 

Since intermediate-range gravitational force has not been discovered in the 

(2) 

solar system, it must be small. The solution of C, has a (1 — 7) factor. The 
intermediate range part of has a (1 — 7)^ factor. From the constraint of 
empirical test (Bertotti et al., 2003), (1 — 7) should be less than 10~^. Hence 
the deviation from Newtonian gravity should be less than 10~^ and A2 is not 
much constrained. Further experiments/observations with better precision will 
be able to measure A2 better. From the metric, we can obtain the geodesic 
equation of motion: 

-T-j = Cat + diPN + d2PN, (77) 
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where ctv is the Newtonian acceleration and it has the form 

3„ = IvTV = -G I 'VK-Pdr' - ^ I Ip^df'. (78) 

2 Jv \r — r'\-^ Q-K Jv \r — r'\-^ 

The full equation of motion will be given and discussed in Dong et al. (2007). 



6 Outlook 



Full 2PN approximation of general scalar-tensor theory of gravity has been 
obtained in a single frame. Multiple-frame studies hke those have been worked 
out for general relativity (Damour Soffel and Xu , 1991, 1992, 1993), for PPN 
formalism with two parameters (Klioner & Soffel, 2000), and for scalar-tensor 
theory (Kopeikin & Vlasov, 2004) in the IPN approximation would be the 
next step to investigate. 2PN approximation of vector-tensor theory of gravity 
(Jacobson & Mattingly , 2001; Luo et al., 2007) and other theories need to be 
worked out to see a more general structure of 2PN approximation. With this 
done, it would be easier to formulate a useful parameterized 2PN formalism 
for testing 2nd order relativistic gravity. 

This research is supported by the National Natural Science Foundation of 
China under Grants No. 10475114 (W.-T. Ni) and No. 10563001 (Y. Xie and 
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